NAisA  leclmical  MeiKorsiidiisii  83694 


310 

KX'^d'-j'-rr 


iY:4e 


k'-T  _ 

JJrUin. -4 <; '} 


P.  L.  N.  Murthy  and  C.  C.  Chsmis 
Lewis  Research  Center 
Cieveiand,  Okio 


DTIC  QUALITY  mSPEU-mu  % 


DYJiAKIC  STRESS  ANALYSIS  OF 
KOTCKED  FlfiER  CCflFOSlTE  FLEXU5AL 
SPECIfiiiis  (K^SA)  ^7  p  iiC  A03/aF  A0 1 

CSCL  11D 


G3/24 


Prepared  for  the 

Seventh  Conference  on  Composite  Materials:  TestiBs-  and  Dei: 

™"“eW-:7  5f  Materia!;  ' 

i  niiaueipLia,  Pennsyh'auia,  April  2-5,  1984 


DEPA';r:-^‘-.-- 

^  i'^sr/cs  ]  ^  ^  '■  -iNSE 

ARRADCOiH, 

07801 


19960226  007 


N84-25770 


Ouclas 

19418 


THIS  DOCUMENT  IS  BEST 
QUALITY  AVAILABLE.  THE 
COPY  FURNISHED  TO  DTIC 
CONTAINED  A  SIGNIFICANT 
NUMBER  OF  PAGES  WHICH  DO 


NOT  REPRODUCE  LEGIBLY. 


I 


TABLE  OF  COfiTEfiTS 

PAGE 


SUmRY  .....  . . 1 

INTRODUCTION  1 

ANALYSIS  . . .....•......!!!!!!!!!!  2 

Specimen  Geometry  ......  .  .  *  ]  2 

Finite  Element  Idealization  . .  2 

Direct  Transient  Response  Analysis  .  3 

Composite  Systems  Analyzed  .  4 

RESULTS  AND  DISCUSSION ..’!!!!  4 

Load  Conditions . 4 

Natural  Frequencies  and  Normal  Hodes  .  4 

Displacement  and  Stress  Wave  Propagation  .  5 

Transient  Response  of  Soiooth  Speci(r;ens .  6 

Static  and  dynamic  Stress  Contours  .  7 

Smooth  Specimen . 7 

Notched  Specitnen  .  . . 7 

GENERAL  DISCUSSION  .  7 

SUNJ1ARY  OF  RESULTS .  g 

REFERENCES  . .  q 


\ 

■ 

I 


It 


P 


1 


E-2152 


DYNAMIC  STRESS  ANALYSIS  OF  SMOOTH  AND  NOTCHED  FIBER 
COMPOSITE  FLEXURAL  SPECIMENS 


P.  L.  N.  Murthy  and  C.  C.  Chamis 
National  Aeronautics  and  Space  Administration 
Lewis  Research  Center 
Cleveland,  Ohio  44135 


SUMMARY 

fiho/  Metalled  analysis  of  the  dynamic  stress  field  In  smooth  and  notched 
fiber  composite  (Charpy-type)  specimens  Is  reported  In  this  paper.  The  analy¬ 
st  Is  performed  with  the  aid  of  ihe  direct  transient  response  analy  Is  tSiu 
tion  sequence  of  MSC/f)ASTRAN.  Three  unidirectional  composites  were  choseS  for 

S-G1ass/Epoxy,  Kevlar/Epoxy  and  T-300/Epoxy  composite 
systems.  The  specimens  are  subjected  to  an  Impact  load  which  Is  modeled  as  a 
triangular  Impulse  with  a  maximum  of  2000  1b  and  a  duration  of  1  ms  The  re- 

rJeak^1oad°T200o1J^  analysis  of  the  specimens  subjected  to 

^atlr  S'  geometry  and  type  of  materials  studied,  the 

static  analysis  results  gave  close  conservative  estimates  for  the  dynamic 

Interesting  Inference  from  the  study  Is  that  the  Impact 
S.G,.ss/Ep„«,  spec„,,„s  sooner  .Pan  .evlar/EpSx, 


IHTROOUCTION 

Smooth  and  notched  flexural  specimen  testing  continues  to  be  popular  In 
the  composites  consrtunlty  for  characterizing  and/or  qualUylnn  fiber  composites 
Reasons  for  the  popularity  are:  (1)  simplicity,  (2)  adaptability  to  SJsl 

P  availability  of  simple  equations  for  data  reduction,  and  (4) 
allablllty  of  written  AST?!  standard  testing  procedures.  In  addition  to 
these  reasons,  flexural  specimen  ter41ng  forces  the  material  t^respofld  like  a 
structure  by  simultaneously  subjecting  It  to  tensile,  compressive  and  shear 
strepes.  Furthermore,  the  notched  flexural  specimen  (Charpy  Impact  test 
renstance!^^^^"^  enables  easy  determination  of  fracture  toughness  and  Impact 

event^'nf^as  ‘Juasl-static  or  dynamic)  Is  a  dynamic 

event  and,  as  such.  Is  a  very  complex  process,  fracture  progression  In  a 

exural  specimen  Is  controlled  by  the  local  dynamic  stress  field  which  Is  also 
complex.  The  local  dynamic  stress  field  Is  chLacterIzed  by  the  stress  wave 
which  are  normal,  shear,  flexural  and  surface  waves.  Each  of  the  dynamic 

de^rL'?"  ®  ^  propagate  this 

descHotlorof^thrin  of  these  phenomena  a  detailed 

description  of  the  local  dynamic  stress  field  Is  a  prerequisite.  This  reoulre-: 

nn*  3  complex  transient  analysis  Involving  direct  time  Integration  as 

standard  modal  synthesis  methods  such  as  those  available  In 
some  general  purpose  finite  element  structural  analysis  codes  Solution  27 

su;"anal;":''1hr:k'?:cti  Is  one  of  the  avanabL“  t1o?s  for 

1  i  '  Th6  obj6Ct1v8  of  tho  prosent  pspor  Is  to  reoort  a 

analysis  of  the  dynamic  stress  field  In  smooth  and  notched  ChaJpy  specimens. 


Although  the  notched  Charpy  test  specimen  has  been  used  for  years  In  test¬ 
ing  metals  and  recently  1n  testing  composites,  no  analysis  has  been  performed 
to  determine  the  detailed  dynamic  stress  state  variation  In  the  notch  vicinity. 
In  general,  the  physical  problem  of  the  notched  Charpy  test  specimen  Is  dynamic 
and  nonlinear;  solution  of  this  problem  Is  difficult.  However,  a  good  first 
order  approximation  may  be  obtained  by  assuming  linear  behavior  and  a  quasi¬ 
static  load.  Reference  1  reports  such  an  analysis  which  uses  the  static  solu¬ 
tion  sequence  of  COSHIC  NASTRAN.  The  main  conclusions  of  reference  1  are;  1) 
the  stress  state  Is  biaxial,  and  2)  the  Charpy  test  specimen  Is  not  suitable 
for  assessing  the  Impact  resistance  of  nonmetalllc  fiber  composites  directly. 

In  the  present  work  the  effort  Is  directed  towards  understanding  the  stress 
wave  propagation  and  the  attendant  dynamic  stress  field  In  fiber  composite 
smooth  and  notched  flexural  specimens. 


ANALYSIS 

In  this  section,  the  specimen  geometry,  the  finite  element  Idealization, 
the  finite  element  analysis  method,  and  the  composite  systems  analyzed  are 
described. 


Specimen  Geometry 

The  geometry  of  the  Charpy  test  specimen  (ASTH  STD  E23-?)  Is  shovjn  In 
figure  1.  As  can  be  seen  In  this  figure,  the  overall  length  of  the  specimen 
Is  2.164  In.  and  the  length  between  supports  Is  1.574  In.  The  specimen  width 
Is  0.394  In.  The  specimen  unnotched  depth  Is  0.394  In.  and  the  depth  at  the 
notch  Is  0.315  In.  The  notch  1s  0.079  In.  deep  and  has  a  45°  opening. 


Finite  Element  Idealization 

Two  finite  element  Idealizations  of  Charpy  test  specimens  --  one  without 
a  notch  (smooth  specimen)  and  one  with  a  notch  were  selected  for  the  present 
study.  The  details  of  the  meshes  are  shown  In  figures  2  and  3.  The  material 
properties  are  assumed  to  be  uniform,  orthotropic,  and  obey  a  linear  stress- 
strain  law  throughout  the  analysis.  In  addition,  the  specimen  Is  assumed  to 
be  In  a  state  of  plane  stress.  The  plane  stress  assumption  Is  also  justified 
from  the  physics  of  the  problem.  The  width  restraints  at  the  notch-tip  are 
negligible  because  of  the  very  low  value  of  the  respective  Poisson's  ratio. 

For  the  present  analysis,  the  plane  stress  assumption  Implies  that  the  stresses 
are  permitted  to  vary  along  the  specimen  length  and  through  the  thickness  but 
not  across  the  width.  This  reduces  the  stresses  to  be  calculated  to  three, 
two  normal  and  one  shear. 

With  these  assumptions,  plane  stress  finite  elements  can  be  used  to  model 
the  Charpy  test  specimen.  For  the  smooth  specimen  all  the  elements  are  quad¬ 
rilateral.  The  notched  specimen  Is  modeled  with  both  triangular  elements  and 
quadrilateral  elements.  The  triangular  elements  are  used  as  transition  ele¬ 
ments  In  the  areas  around  the  the  supports,  the  load  application  point  and  the 
notch.  These  are  the  regions  where  maximum  stress  concentrations  are  expected 
to  occur  and  therefore  are  provided  with  a  finer  mesh.  The  boundary  conditions 
prescribed  are  such  that  the  node  at  the  left  support  Is  constrained  from  x,  y 
and  z  displacements,  and  the  node  at  the  right  support  Is  constrained  from  y 


and  z  displacements.  In  addition,  for  the  notched  specimen,  three  nodes  clos¬ 
est  to  the  right  support  are  constrained  from  displacement  In  the  y  direction. 
The  specimens  are  subjected  to  an  Impulse  loading.  The  form  of  Impact  Is  a 
triangular  function  with  peak  load  of  2000  lb  occuring  at  500  ys  of  a  total 
contact  time  of  1000  ys. 

The  statistics  of  the  finite  element  representation  for  the  smooth  and 
notched  specimens  are  as  follows: 


1)  Smooth  Specimens. 

Number  of  nodes  or  grid  points  1147 

Number  of  displacement  degrees  of  freedom  (OOF)  (2  degrees 
of  freedom  per  node)  2294 

Number  of  quadrilateral  plate  elements  (CQUAD4)  1080 

DOF  eliminated  using  the  boundary  conditions  (2  from  u  »  0 
and  V  =  0  at  the  left  support  and  1  from  v=0  at  the 

right  support)  3 

Number  of  free  OOF  (2294  -  3)  2291 

2)  Notched  Specimens. 

Number  of  nodes  or  grid  points  056 

Number  of  displacement  degrees  of  freedom  (DOF)  1312 

Number  of  quadrilateral  plate  elements  (CQUA04)  544 

Number  of  triangular  plate  elements  (CTRIA3)  82 

Total  number  of  elements  (544  *■  82)  626 

DOF  eliminated  using  the  boundary  conditions  (2  from  u  =  0 
and  V  =  0  at  the  left  support  and  3  from  v  ^  0  for  three 
nodes  at  the  right  support)  5 

Number  of  free  DOF  (1312  -  5)  1307 


Finite  Element  Analysis  Method 

The  MSC/NASTRAN  general  purpose  structural  analysis  finite  element  com¬ 
puter  program  Is  used  for  the  finite  element  analysis.  The  specific  elements 
used  are  Identified  as  CTRIA3  and  C0UAD4.  They  are  Isoparametric  constant 
strain  elements.  NASTRAN  obtains  the  solution  using  a  displacement  formulation 
via  rigid  format  solution  sequence  No.  27.  This  solution  sequence  employs  a 
direct  time  Integration  scheme  to  obtain  the  transient  response  of  a  structure. 

The  solution  sequence  No.  27  of  MSC/NASTRAN  uses  the  Integration  algorithm 
based  upon  the  Newmark  Beta  method  (ref.  2).  It  provides  stable  results  for 
the  widest  possible  spectrum  of  practical  problems  without  sacrificing  either 
accuracy  or  efficiency.  For  complete  details  of  the  MSC/NASTRAN  analysis 
(ref.  3)  should  be  consulted.  A  brief  description  Is  given  In  the  following 
paragraphs.  j 

The  differential  equations  of  a  linear  structural  problem  may  be  written 
In  the  general  matrix  form 

[p^M  f  pB  ♦  K]  (u)  =  (F)  (1) 

where  p  =  d/dt.  H  and  K  are  the  mass  and  stiffness  matrices,  B  Is  asso¬ 
ciated  with  the  damping  matrix,  u  Is  the  vector  of  displacements  and  F  Is 
the  load  vector.  The  numerical  Integration  Is  achieved  by  replacing  p2  and 
p  by  finite  difference  operators  and  then  using  explicit  Integration.  The 
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outputs  of  the  transient  analyst?  module  Include  velocities,  and  accelerations 
as  well  as  displacements.  The  output  can  be  requested  at  even  multiples  of 
the  Integration  time  step.  This  feature  affords  some  economy  In  output  data 
preparation  In  cases  where  small  time  step  Is  needed  for  greater  accuracy. 


Composite  Systems  Analyzed 

Flexural  (Charpy-type)  test  specimens  made  from  three  typical  composite 
systems  are  analyzed.  They  are:  T-300/Epoxy,  Kevlar/Epoxy,  and  S-Glass/Epoxy 
composites.  The  specimens  are  all  unidirectional  composites  with  the  fibers 
parallel  to  the  length  (x-axIs,  fig.  1)  of  the  specimen. 

The  plane  stress-strain  relationship  (stiffness)  coefficients  required  to 
Input  to  NASTRAN  are  surmarlzed  In  table  I.  These  properties  are  obtained  by 
using  the  resident  data-bank  In  the  composite  mechanics  computer  code  ICAN 
(ref.  4).  The  relationships  between  the  NASTRAN  stiffness  coefficients  (G's), 
and  the  usual  engineering  constants  are: 


«11 

=  Si/(1  - 

*'112  ’'121^ 

(2) 

®12 

=  ''121  Si 

=  *'112  S2  °  Si 

(3) 

G22 

=  E, 22/(1  - 

*’112  *'121^ 

(4) 

S3 

"  *"112 

(5) 

S3 

""  S3  Si 

=  S2  =  ° 

(6) 

The  notation  In  equations  (2)  to  (5)  Is  as  follows:  E^n  denotes  the  longi¬ 
tudinal  modulus,  Ei22  transverse  modulus,  Gjn2  the  shear  modulus,  ui']2 
the  major  Poisson's  and  »j^21  minor  Poisson's  ratios.  For  an  elastic 
material  the  two  Poisson's  ratios  are  related  by  the  well  known  relation 

”il21  "  ’■'!n2  ^1122'^^iLll 
RESULTS  AND  DISCUSSION 
Load  conditions 

Two  types  of  load  conditions  are  used  In  obtaining  the  results.  The  first 
Is  a  static  loading  where  a  force  of  2000  1b  Is  applied  at  the  center  of  the 
specimen  on  the  top  surface.  The  second  Is  an  Impulse  loading.  This  Is 
modeled  as  a  triangular  pulse  with  a  peak  value  of  2000  lb  In  the  middle.  The 
pulse  Is  modeled  to  last  for  1000  ps.  The  transient  response  Is,  however, 
obtained  for  three  contact  time  periods  (I.e.,  3  ms).  A  separate  normal  modes 
analysis  Is  used  to  determine  the  first  five  natural  frequencies. 


Natural  Frequencies  and  Normal  Modes 

The  NASTRAN  normal  modes  analysis  module  (Solution  3)  1s  utilized  to  de¬ 
termine  the  first  five  natural  frequencies  and  the  associated  mode  shapes  for 
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the  Ch-arpy  test  specimens.  They  are  shown  In  figure  4.  The  results  are  shown 
•n  table  II  for  S-61ass/Epoxy  specimens.  The  frequencies  are  used  to  determine 
the  time  periods  which  aid  In  determining  the  Integration  time  step  for  the 
transient  analysis.  A  time  step  of  5  ps  Is  chosen  for  the  transient  analysis 
based  upon  the  time  periods  shown  In  table  II. 


Displacement  and  Stress  Wave  Propagation 

The  bulk  wave  and  shear  wave  velocities  are  normally  much  higher  compared 
to  the  flexural  wave  velocities.  In  order  to  capture  the  characteristics  of 
propagation  of  these  waves  a  much  smaller  time  step  of  Integration  (0.1  ps)  Is 
chosen.  The  output  of  displacements  and  stress  contours  at  various  time  Inter¬ 
vals  are  saved  and  displayed  graphically  In  figures  5  to  10.  Figures  5  to  7 
depict  the  dynamic  displacement  propagation  In  S-Glass/Epoxy,  Kevlar/Epoxy  and 
T-300/Epoxy  specimens.  Two  bulk  wave  velocity  parameters  and  one  shear  wave 
velocity  parameter  are  defined  below  to  aid  the  following  discussion  of  the 
displacement  wave  propagation  results: 


Bll  " 

6n/P 

(8) 

B22  "" 

G22/P 

(9) 

SI  2  * 

G12/P 

(10) 

where  p  Is  the  mass  density  of  the  material. 

The  computed  values  of  Cg-)-!,  C322>  *-$12  shown  In  table  III  for 

the  three  composite  systems  under  study.  The  velocities  are  expressed  In 
In./pS.  The  trend  indicated  by  Cq22  ^^r  the  three  materials  (the  transverse 
shock  wave  travels  fastest  In  S-Glass/Epoxy  and  slowest  In  Kevlar/Epoxy)  Is 
seen  clearly  In  figures  5  to  7  specifically  the  frames  after  3  and  5  pS.  A 
rough  estimate  of  the  normal  wave  velocity  can  be  obtained  by  counting  the 
number  of  elements  that  appear  to  be  affected  by  the  Impact  from  the  figures  5 
to  7.  The  normal  wave  velocity  estimates  from  the  1  and  3  sec  frames  are  shown 
below: 


Composite  system 

Number  of 

1  psec 

elements 

3  psec 

Velocity 

1  psec  3  psec 

Average 

S-Glass/Epoxy 

10 

25 

0.1313 

0.1094 

0.1204 

T-300/£poxy 

7 

16 

.0788 

.0700 

.0744 

Kevlar/Epoxy 

5 

12 

.0657 

.0520 

.0589 

These  values  are  In  close  agreement  with  the  theoretical  values  shown  In 
table  III  under  Cgp?-  The  same  trend  Is  also  seen  for  the  3  and  5  p  sec 
frames.  It  can  be  concluded  that  the  Initial  shock  travels  with  the  bulk  wave 
velocity  Cg22  alnng  the  direction  of  Impact. 

Once  the  normal  shock  reaches  the  bottom  of  the  beam,  the  wavefront  ap¬ 
pears  to  be  moving  1n  the  longitudinal  direcclon  forming  a  flexural  wave.  The 
velocities  of  the  waves  traveling  In  the  longitudinal  direction  can  also  be 
determined  approximately  with  the  same  technique  mentioned  earlier.  The  fol¬ 
lowing  are  the  details  for  the  frame  after  13  ps: 


Compos Vte  system  Number  &f 
elements 


Veloct-ty 


S-Glass/Epoxy 

T-300/£poxy 

Kevlar/Epoxy 


11  0.0510  In/psec 

9  .0416  In/psec 

7  .0320  tn/jjsec 


The  above  velocities  appear  to  have  the  same  trend  as  depicted  by  the  wave 
velocity  parameter  Cs'j2  shown  1n  table  III.  However,  the  waves  along  the 
longitudinal  direction  appear  to  move  significantly  slower  than  that  Indicated 
by  Csi2.  This  Is  probably  due  to  the  coupling  between  the  flexural  wave  and 
the  shear  wave.  The  flexural  wave  velocity  Is  significantly  slower  than  the 
shear  wave  velocity.  For  example,  the  smooth  S-Glass.^Epoxy  specimen  under 
study  has  a  flexural  wave  velocity  given  by 


Cp  =  2lf  =  0.0137  In/jisec 


where  i.  Is  the  length  between  the  supports,  f  1s  the  first  fundamental 
frequency.  (It  Is  assumed  that  the  wave  number  Is  1  and  the  beam  deflects  Into 
a  half  wave.) 

The  transient  stress  response  after  one  percent  (10  ps)  of  the  contact 
time  are  shown  In  figures  8  to  10  In  the  form  of  stress  contours.  All  stresses 
are  localized  at  this  early  time.  The  longitudinal  stress  (<Jn)  in  the 
Kevlar/Epoxy  and  T-300/Epoxy  composites  Is  of  about  the  same  magnitude  while 
that  In  the  S-Glass/Epoxy  Is  about  half  as  much.  The  normal  stress  (022) 
shear  stress  (0^2)  are  of  about  the  same  magnitude  for  all  three  composite  sys¬ 
tems.  It  appears  from  these  stress  results  that,  under  the  same  Impact  condi¬ 
tions.  the  stress  In  the  Kevlar/Epoxy  and  the  T-300/Epoxy  will  reach  fiber 
fracture  stress  magnitudes  considerably  earlier  (about  half  the  time)  than  In 
the  S- Glass/Epoxy  composite.  Two  Implications  follow  relative  to  the  same 
stress  magnitude:  (1)  the  S-Glass/Epoxy  composite  will  sustain  greater  Impact 
load  prior  to  fracture  than  the  T-300/Epoxy,  and  (2)  the  rapid  compressive 
stress  built  up  will  cause  longitudinal  compression  failure  accompanied  by 
substantial  bending  deflection  In  the  Kevlar/Epoxy  composite  thus  Increasing 
the  Impact  required  to  Induce  fracture. 


Transient  Response  of  Smooth  Specimens 

The  displacements,  velocities,  accelerations  and  stresses  are  obtained  for 
a  total  contact  time  of  1000  iis  for  S-Glass/Epoxy  smooth  specimens.  In  these 
computations  a  time  step  of  5  \i$  Is  used.  The  output  Is  saved  for  every  two 
time  steps.  Two  points  A  and  B  as  shown  In  figures  2  and  3  are  selected  for 
study.  A  Is  the  load  point  and  B  1^  the  opposite  point  at  the  bottom  for 
smooth  and  at  the  notch-tip  for  notched  specimens. 

The  transient  response  results  appear  In  figures  11  to  16.  Figures  11  to 
13  shew  the  longitudinal  (axial)  and  the  transverse  (bending/flexural)  com¬ 
ponents  of  displacement,  velocity  and  acceleration  of  point  A  plotted  against 
time.  The  corresponding  results  for  point  B  are  not  shown  as  they  are  similar 
to  point  A  response  both  qualitatively  and  quantitatively.  It  can  be  concluded 
from  these  figures  that  the  response  Is  primarily  In  the  first  flexural  mode. 
For  example  (from  fig.  11)  the  number  of  cycles  In  1  ms  Is  counted  as  a  little 
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over  6.  From  table  II  the  first  fundamental  time  period  Is  161.3  ps  whlc'h 
Implies  that  6.2  (1000/161.3)  cycles  of  response  In  first  mode  can  be  expected. 

The  stress  response  results  are  shown  graphically  In  figures  11  to  16. 

Each  figure  contains  the  stress  response  of  elements  near  points  A  and  B.  The 
longitudinal  stress  behavior  appear  to  be  primarily  the  first  flexural  mode 
response.  However,  the  magnitude  of  stress  near  point  A  Is  about  one  and  a 
half  times  greater  than  the  magnitude  of  stress  near  point  B.  This  Is  to  be 
expected  because  of  stress  concentration  near  the  load  point  A.  The  transverse 
normal  (022)  ^nd  the  shear  (oip)  stress  Increase  linearly  to  a  peak  when  the 
load  takes  the  maximum  value  of  2000  1b  and  then  decrease  linearly  to  a  zero. 
The  corresponding  stresses  for  the  element  near  point  B  are  Insignificant  and 
therefore  can  not  be  shown  distinctly  when  drawn  to  the  same  scale. 

The  response  after  the  load  removal  Is  not  shown  In  the  above  figures. 

The  computations,  however,  are  conducted  for  three  contact  time  periods.  It 
Is  found  that,  for  all  practical  purposes,  the  specimen  remains  static.  A 
typical  response  curve  Is  shown  In  figure  17  for  point  A. 


Static  and  Dynamic  Stress  Contours 

Smooth  specimen.  -  The  static  stress  contours  under  a  load  of  2000  1b  and 
the  dynamic  stress  contours  at  the  peak  load  of  2000  lb  are  shown  In  figures 
18  to  20.  They  appear  to  be  Identical.  The  peak  values  for  the  dynamic  load 
case,  however,  are  slightly  lower  than  those  for  the  static  case.  While  the 
maxima  for  the  longitudinal  stress  (o]^}  occur  at  the  center  of  top  and  the 
bottom,  the  transverse  normal  stress  and  the  transverse  shear  stress  have  the 
maxima  almost  adjacent  to  each  other  near  the  load  application  point.  The 
transverse  shear  stress  and  the  axial  stress  distribution  Is  seen  to  approach 
the  classical  Euler-Bernoull 1  theory  predictions  away  from  the  load  application 
point.  Steep  stress  gradients  are  observed  near  the  load  application  point. 

Notched  specimen.  -  The  results  for  notched  S-Glass/Epoxy  specimen  are 
shown  In  figures  21  to  23.  The  static  and  dynamic  stress  contours  at  the  peak 
load  show  similar  trend  as  observed  1n  the  case  of  smooth  specimens.  Figure  24 
shows  the  stress  Intensity  near  notch-tip  In  relation  to  the  far  field 
stresses.  This  distance  from  the  notch-tip  Is  measured  along  the  longitudinal 
direction  towards  the  left  support.  It  Is  seen  that  all  the  three  stresses 
attain  very  high  peaks.  Indicating  severe  local  stress  Intensities  near  the 
notch-tip.  As  the  stress  allowables  for  022  °12  generally  an  order 

of  magnitude  lower  than  that  for  o^-j,  one  can  expect  a  matrix  Initiated  fail¬ 
ure  followed  by  fiber  fractures  at  this  location. 


GENERAL  DISCUSSION 

The  response  analysis  of  smooth  and  notched  Charpy-type  flexural  specimens 
subjected  to  a  triangular  Impulsive  load  provides  insight  Into  the  nature  of 
stress  wave  propagation,  failure  mechanisms  and  the  relation  between  the  static 
and  dynamic  responses.  The  contact  time  of  the  Impact  load  Is  approximately  6 
time  periods  of  the  first  flexural  mode.  The  bulk  wave  velocities  and 

Cb22»  the  shear  j>/ave  velocity  Cs')2  very  high  compared  to  the  flexural 
wave  velocity  for  the  materials  and  the  geometry  under  study.  Hence,  the  re¬ 
sponse  Is  observed  at  early  times  of  the  Impact  event  (of  the  order  of  a  micro¬ 
second).  It  shows  that  the  Initial  waves  travel  with  a  bulk  wave  velocity 
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Cb22  reach  the  bottom  surface  of  the  specimen.  The  waves  then 

appear  to  travel  longitudinally  towards  the  supports.  The  approximate  calcu¬ 
lations  based  upon  the  deformations  at  various  times  Indicated  that  the  wave 
velocity  Is  significantly  lower  than  the  shear  wave  velocity  Csi2  but  much 
higher  than  the  flexural  wave  velocity.  Ihe  wavefront  appears  to  Induce  a 
flexural  wave  as  It  progresses  towards  the  supports. 

At  longer  times,  the  response  Is  primarily  In  the  first  flexural  mode.  A 
comparison  with  the  static  response  Indicates  that  little  or  no  difference 
exists  In  the  magnitude  for  stresses  at  the  peak  load.  The  static  predictions, 
however,  are  on  the  conservative  side. 

The  response  of  notched  fiber  composite  specimens  show  severe  stress  In¬ 
tensities  near  the  notch  region  Indicating  local  failures  In  shear  and  trans¬ 
verse  tension  Initially  followed  by  fiber  fractures  since  the  longitudinal 
stresses  approach  fracture  stress  magnitudes. 

The  response  after  the  load  removal  reduces  to  mere  noise;  the  .structure 
remains  practically  static. 


SUMMARY  OF  RESULTS 

The  results  of  the  study  on  static,  dynamic  and  transient  response  of 
smooth  and  notched,  fiber  composite,  Charpy  specimens  are  listed  below: 

(1)  The  wave  propagation  velocities  can  be  estimated  from  the  early  time 
displacement  propagation  response;  the  estimates  are  In  fair  agreement  with 
the  theoretical  predictions. 

(2)  The  first  five  flexural  mode  shapes  and  frequencies  show  coupling 
effects  from  the  thickness  stretch  and  the  transverse  shear  type  modes;  the 
frequencies  are  not  Integer  multiples  of  the  fundamental  frequency. 

(3)  The  dynamic  and  static  peak  load  stress  contours  are  almost  Identi¬ 
cal.  The  static  peak  load  stress  magnitudes  are  slightly  higher.  Stress  pre¬ 
dictions  based  upon  a  quasi-static  approach  lead  to  conservative  estimates. 

(4)  The  transient  displacement  and  velocity  response  appear  to  be  pri¬ 
marily  In  the  first  flexural  mode.  The  acceleration  response  shows  contribu¬ 
tions  from  higher  modes. 

(5)  The  structure  responds  for  the  load  duration  time  period  only.  The 
steady  state  response  after  the  load  removal  appears  to  be  negligible. 

(6)  The  transverse  normal  and  shear  stresses  vary  linearly  with  time  and 
follow  the  load  path. 

(7)  The  notch-tip  region  develops  severe  stress  concentrations  and  any 
of  the  three  stresses  could  cause  or  Initiate  a  failure. 

(8)  Based  upon  the  transient  stress  response.  It  appears  that  the  fall- 
•irifs  1n  Charpy  specimens  are  Initiated  at  the  notch-tip  by  the  shear  stresses 
and  the  transverse  tensile  stresses  followed  by  fiber  fractures.  However,  the 
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failures  in  smooth  specimens  are  probably  inltlatpd  by  a  combined  stress  state 
near  the  load  application  point  followed  by  local  Interply  delaminations. 
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ORICSMAL 
OF  POOR  QUALITY 


TABLE  I.  .  NASTRAN  PLANE  STRESS  ~  STRAIN  RELATIONSHIPS 


Composite 

Mass  density 

Stress-strain  coefficient 

Orthotropy  Ra 

tios 

system 

(xIO*) 

(1b.sec^/1n*) 

*^12  *  *^2) 

®22 

(xIO^  psi) 

Sa 

''ll/''22 

S-Glass 

/Epoxy 

1.813 

8.7900 

1.1700 

3.2600 

1.2700 

2.70 

6.92 

Kevlar 

/Epoxy 

1.249 

12.385 

.20082 

.57377 

.31635 

21.59 

39.15 

T-300 

/Epoxy 

1.405  1 

17.896 

j 

.30275 

1.1318 

i 

.50371 

15.81 

35.18 

TABLE  II.  -  KATURAL  FREOUEKCIES  AND  TIHE 
PERIODS  FOR  S-6LASS/EP0XY  5HOOTH  AND 
NOTCHED  SPECIMENS 


Specimen 

type 

Mode 

number 

Frequency 

(cycles/sec) 

Tlfne  period 
(y$) 

Smooth 

1 

6200.67 

161.3 

2 

11659.54 

85.8 

3 

19993.02 

50.2 

4 

31062.8 

32.2 

5 

43313.02 

22.8 

Notched 

1 

9030.16  1 

110.7 

2 

14631.33  ! 

68.3 

3 

24596.56  1 

40.7 

4 

27775.86  1 

1  36.0 

5 

36685.27 

27.3 

TABLE  III.  -  HAVE  VELOCITY  PARAMETERS 
FOR  THE  SELECTED  COMPOSITE  SYSTEMS 


jr  e  ter 

Mater 

Sn 

^B22 

St2 

S-Glass/Cpoxy 

0.2201 

0.1341 

0.0837 

Kevlar/Epoxy 

.3149  1 

.0678  1 

.0503  1 

T-300/epox> 

.3569 

.0898  1 
1 

.0602  i 
_ i 

ORICSNAt  PA»^S 

OF  POO^.  Q'J  ALTY 


L  f  •  6200. 67  Hz.  1.  f .  9030. 157  Hz. 


4.  f-  31062.80  Hz.  4.  f  •  27775.  ^  Hz 


5.  f »  43913. 02  Hz.  5.  f .  35585  27  ^ 


Figure  4.  -  Natural  frequencies  and  mode  shapes  for 
S-giass/epoxy  smooth  and  notched  flexural  specimens. 
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Figure  5.  -  Dynamic  dispUcemsnt  propagation  in  S-glass/ 
epoxy  composite  flexural  specimen,  (Triangular  impulse: 
2000  lb  roaxim  jm  1000  usee  duration. ) 
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Figure  6.  -  Dynamic  displacement  propagation  in  Kevlar/ 
epoxy  composite  flexural  specimen.  (Triangujar  impulse 
2000  lb  maximum  1000  usee  duration. ) 
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Figure  7.  •  Dynamic  displacement  propagation  in  T~300/ 
epoxy  eomposite  flexural  specimen.  (Triangular  impulse 
2(XX)  lb  maximum  1000  usee  duration. ) 
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SVMSOL  VALUES,  p$i 

1  -1.328S99E*03 

2  -1.  143674E+03 

3  -9. 5844g9f +02 
'  4  -7.732237E+02 

5  -5.S79985eK^ 

6  -A027733£+02 

7  '2AVymm 
B  -3.23229EE+01 
9  L  529022L4-02 

10  3.3S1274£f02 


''  L 


9__- 

_  10  . 

to)  aj2  Contcxirs. 


1  -L052444E+O3 

2  -9. 068754E+02 

3  -7.613066E+02 

4  -9. !  57379c +02 
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6  3.246003E+02 

7  L7S^3i6E+02 

8  3. 346?a3E+01 

9  l.l2l05->£ +02 
10  Z  576747L+02 


1  -5. 760!Sa(?£+02 

2  -A  970201 E+02 

3  -4.120S02L+02 

4  -3.689U4E+02 

5  -^998825E+0^ 

6  -1, 808336E+02 

7  -1.017847E+02 

8  *2. 273584E+01 

9  5,631304£+01 

10  L  353619E+02 


(c)  ©22  Contours, 


Figure  a  -  Dynamic  stress  contours  in 
s-gass/epoty  composite  flexural 
specimen  after  1  percent  contact  tiraa 


ORIGINAL  PAGE  fS 
OF  POOR  QUALITY 


SYMBOL  VALUES,  psi 

1  '^  362945E+03 

2  -ao&5930E+03 

3  'L768912E+03 

4  -1,47I895E+03 

5  -L 174377E+03 

6  *a77fi5e9E+02 

7  *5.30S41lE+02 

8  “^838?33£+0^ 

9  1 31M53£:+0I 
10  3. 102323E+02 


1  '1.9I3533E+03 

2  -a  7752i6E+02 

3  “7. 41 5100E+02 

4  -6.054%4£+02 

5  -4694^2t+02 

6  -3.334752E+0? 

7  -1. 97^35f+02 

8  145193E+01 

9  7. 45596^+01 
10  2. 105713f+02 


1  7.77199IE+02 

2  94&C?70f+02 

3  -i920i60£+02 

4  -4. 994229E+02 

5  '<065309£+02 

6  -3.  l423&aE+02 

7  *^216468E+02 

8  -1,290S4^^02 

9  '•3L646273E+C1 
10  5,612931E+01 


after  1  percent  of  the  r antact  time. 
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?  6  -7.38722U^ 

!  7  <  46430^*02 
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9  1.  381526E+02 
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(s)  Oj^  Contours. 


1  'l,019U7r+03 

2  -&  8168188+02 
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6  -3.  3194268+02 

7  -1.9450788+02 
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9  a038^1Si>:+01 
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1  -7.0269858+02 

2  -6.  1789508+02 

3  8309148+02 

4  -4  4828788+02 

5  -3.6348428+02 

6  -2  7863068+02 

7  -1.9387718  +02 

8  -1.0907258+02 

9  -2. 4269918+01 
10  6.0533678  +01 


Figure  10.  -  Dynamic  stress  I'^ontours  in 
T-300/epo>:y  composite  flexural  specimen 
after  1  percent  t  -ntact  time. 


TRANSVERSE  OISPLACEMlNT.  V.  in.  LONGITUDINAL  DISPLACEMENT,  u,  in 
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Specimen. 


Hgure  U  -  Axial  stress  re¬ 
sponse  of  top  and  bottom  elements 
near  points  A  and  B  for  S-glass/ 
epoxy  smooth  flexural  specimen. 


TRANSVERSE  NORAWL  STRESS, 


ORIGINAL  PAGE 
OF  POOR  QUAimf 
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FJgure  15,  -Transverse normal 
stress  Cyy  response  of  top  and 
txjttom  elements  near  points  A 
‘  and  B  for  S-^fass/epcsy  smooth 
flexural  specimen. 


TliV£,  msec 


Figure  16.  -  Transverse  shear  stress 

response  of  top  and  bottom 
elements  near  points  A  and  8  for 
S-glass/epDxy  smooth  flexural 
speciRj'u'H, 


,TIM£,  msec 

Figure  17.  -  Displacement  response 
of  fwint  A  for  three  con^ct  time 
periods  of  S-^lass/epoxy  smooth 
flexural  specimen. 
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(a!  Static  stress  contours. 
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tat  Static  stress  contoyrs, 


1  -1. 4&3251EtG5 

2  -1. 205350E+05 

3  -9. 3044S4Etruj 

4  -6.  555471  £+04 

5  -3. 806458E+04 

6  -L057445E-H34 

7  1.691569E+04 

8  4  4«)582£+04 

9  7. 189595E+04 
10  9. 93a6CS£+(>4 


tbl  Dynam'C  stress  contours. 


Figure  la  -  Coroparison  of  static  an-s  dynamic 
axial  stress  contours  at  pee:<  i.aid  (or 
“9  ass/epuxy  s-nuith  flexural  specimens. 
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5  -Z  3503332+04 
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ngure  .g.  -  Comparison  of  static  and  dynamic 
transvers?  normal  stress  o...  contours  at 
PMk  load  for  S-olass/epoxy  Sfr.ootn  flexural 
specimen. 
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<a)  Static  stress  ccifours. 
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(b)  Dynamic  stress  contours. 


Figure  20.  -  Comparison  of  static  and  dynamic 
transverse  shear  stress  contours  at 
peak  load  for  S-glass/epoxy  smooth  fiexural 
specimen. 
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<ai  Static  stress  corttours. 
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2  Z 137141E+C4 

3  6.359>a2r*04 

4  L.06&1«E+05 

5  L  492362E+05 

6  L918579E+05 

7  Z  34479X+05 

8  2.771011E+05 

9  3. 197227E+05 
10  3.623443E+05 


(b)  Dynamic  stress  contours. 


Figure  21,  r  Comparison  of  static  and 
dynamic  axial  stress  contours 
near  the  notch  region  at  psk  load 
for  S-glass/epaxy  notched  flexural 
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(a)  Slatic  stress  contours. 
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(b)  Dyna.Tiic  stress  contours. 


Figure  22.  -  Comparison  of  static  and 
dynamic  transverse  normal  stress 
Oy^  contours  near  the  notch  region 
at^ak  load  for  S-glass/epory 
notched  flexural  specimen. 
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(a)  Static  stress  contours. 
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(b)  Dynamic  stress  contours. 


Figure  23.  -  Comparison  of  static  and 
dynamic  transverse  shear  stress 
0^  contours  near  the  notch  region 
af^k  load  for  S-glass/epoxy 
notched  flexural  specimen. 


